航天器运输包装箱仿真验证平台技术研究

针对数量有限的物理跑车试验无法满足减振与保温性能测试需求的问题，提出一套航天器运输包装箱动力学与热学仿真验证方法，包括：建立适用于包装箱系统的刚柔耦合多体动力学系统，通过结合线路条件测试生成的动力学系统外部激励，实现减振性能虚拟跑车测试；建立基于计算流体力学的包装箱热学模型，通过模拟自然对流和空调控制，实现包装箱保温性能虚拟跑车测试；基于C/S架构和导航式流程设计思想，建立航天器运输包装箱仿真验证平台，通过实际案例证明该平台仿真结果与实际跑车测试数据具有较高的一致性。     

内装式空射运载火箭重力出舱机箭耦合动力学分析

研究了以运输机为平台的内装式运载火箭空射过程载机和火箭的耦合动力学建模。建模针对两个阶段：第一阶段，火箭固定于载机机舱内，两者构成一个整体，按照普通刚体的力学方法处理；第二阶段，火箭解锁后，沿着舱内的发射筒向外滑行，载机和火箭形成两刚体相互作用的耦合系统，基于牛顿-欧拉法建立系统动力学模型。载机在空射火箭过程中，油门和升降舵满偏，在加速前飞的同时拉大姿态俯仰角，火箭在自身重力分量和惯性力的作用下，沿着机舱内的发射筒加速向外滑行，直至与载机分离。数值仿真分析了空射过程载机的重要力学参数的变化过程，验证了载机操控策略的可行性和安全性，可为未来中国空射运载火箭技术设计提供数据参考。     

Effect of bristle pack position on the rotordynamic characteristics of brush-labyrinth seals at various operating conditions

Over the last few decades, the research on the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals is not sufficient. To this end, two kinds of brush-labyrinth seals for the bristle pack element installed upstream of the labyrinth teeth named BSU and installed downstream of the labyrinth teeth called BSD were used to investigate the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals. Using the numerical model combining the porous medium model and the whirling rotor method, the rotordynamic characteristics of the BSU and BSD at various operating conditions including four kinds of pressure ratios, five kinds of inlet preswirl speeds and four kinds of rotor spinning speeds were conducted. The obtained results show that the effects of operating conditions on rotordynamic coefficients for the different seal configurations are different. The direct stiffness, cross-coupled stiffness and direct damping of the BSU are lower than those of the BSD. The rotordynamic coefficients of the BSU are more insensitive to the operating conditions variation. From the perspective of the seal stability, the BSU is a better brush-labyrinth seal configuration at high pressure ratio, high positive preswirl or high rotor spinning speed conditions. While in the case of low pressure ratio, low positive preswirl or low rotor spinning speed conditions, the BSD is a better choice.     

Multi-faults diagnosis of rolling bearings via adaptive customization of flexible analytical wavelet bases

Multi-faults detection is a challenge for rolling bearings due to the mode mixture and coupling of multiple fault features, as well as its easy burying in the complex, non-stationary structural vibrations and strong background noises. In this paper, a method based on the flexible analytical wavelet transform (FAWT) possessing fractional scaling and translation factors is proposed to identify multiple faults occurred in different components of rolling bearings. During the route of the proposed method, the proper FAWT bases are constructed via genetic optimization algorithm (GA) based on maximizing the spectral correlated kurtosis (SCK) which is firstly presented and proved to be efficient and effective in indicating interested fault mode. Via using the customized FAWT bases for each interested fault mode, the original vibration measurements are decomposed into fine frequency subbands, and the sensitive subband which enhances the signal-to-noise ratio (SNR) is selected to exhibit the fault signature on its envelope spectrum. The proposed method is tested via simulated signals, and applied to analyze the experimental vibration measurements from the running roller bearings subjected to outrace, inner-race and roller defects. The analysis results validate the effectiveness of the proposed method in identifying multi-faults occurred in different components of rolling bearings.     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

Development and performance evaluation of tropospheric scintillation model on Ku-band satellite link over Akure,Nigeria

Existing amplitude scintillation prediction models often perform less satisfactorily when deployed outside the regions where they were formulated. This necessitates the need to evaluate the performance of scintillation models developed in one region using data data from other regions while documenting their relative errors. Due to its variation with elevation angle, frequency, other link parameters and meteorological factors, we employed three years (January 2016 to December 2018) of concurrently measured satellite radio beacons and tropospheric weather parameters to develop a location-based amplitude scintillation prediction model over the Earth-space path of Akure (7.17^oN, 5.18^oE), South-western Nigeria. The satellite beacon measurement used Tektronix Y400 NetTek Analyzer at 1 s integration time while meteorological parameters, namely; temperature, pressure and relative humidity were measured using Davis Vantage Vue weather station at 1 min integration time. Comparative study of the model’s performance with nine (9) existing scintillation prediction models indicates that the best and worst performing models, in terms of root mean square error (RMSE), are the Statistical Temperature and Refractivity (STN) and direct physical and statistical prediction (DPSP) models with values 11.48 and 51.03 respectively. Also, worst month analysis indicates that April, with respective enhancement and fade values of 0.88 and 0.90 dB for 0.01% exceedance, is the overall worst calendar month for amplitude scintillation.     

机载托架自动校准系统的设计与实现

为了解决机载托架传统校准方法精度低、耗时、费力等问题,设计了一种高精度、高效率并且操作简单的成品托架自动校准系统。该系统运用调平精度较高的"循环多次"最高点不动调平方法,建立高精度机载托架校准的静力学数学模型,运用VC与Matlab的COM接口编程技术开发机载托架自动校准系统,并建立良好的人机交互界面。试验结果表明,利用该方法对机载托架进行校准,无论是校准精度还是校准效率都较传统方法有很大提高。     

多梁复合结构的动力学建模及全局解析振型分析

多梁复合结构广泛用于工程领域，其动力学特性对相关装备的工作状况有着重要影响，建立准确适用的动力学模型对后续的动力学响应参数影响规律研究有决定性作用。针对工程中常见的多梁复合结构，给出了位移场的一般表达式，利用哈密顿原理得出多梁复合结构的运动控制方程和边界条件，对其全局解析振型进行理论计算，给出了计算多梁复合结构数值振型的一般过程，基于数值分析和实验分析以L型梁和Z型梁为例利用频率相对误差、振型定性比较和MAC定量分析说明了全局解析振型的求解验证方法，利用伽辽金方法和全局解析振型给出了截断原系统得到离散系统运动方程的一般形式。提出的多梁复合结构的建模方法和全局解析振型求解验证方法，对工程中此类结构动力学响应的参数影响规律研究有重要意义。     

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符；反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。